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Abstract. Electronic excitations in a graphitic monolayer (graphene) in the long- 
wavelength approximation are characterized by the linear dispersion law, representing 
a unique example of the really two-dimensional " ultrarelativistic" fermionic system 
which in the presence of topological defects possesses rather unusual properties. A 
disclination that rolls up a graphitic sheet into a nanocone is described by a pointlike 
pseudomagnetic vortex at the apex of the cone, and the flux of the vortex is related to 
the deficit angle of the conical surface. A general theory of planar relativistic fermionic 
systems in the singular vortex background is employed, and we derive the analytical 
expression for the charge which is induced at finite temperature on some graphitic 
nanocones. 
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1. Introduction 

Carbon nanostructures are intensively studied both experimentally and theoretically, in 
view of their expected important implications for the development of electronic devices, 
flat panel displays, nanoswitches, etc. (see, e.g. [T]). The recent synthesis of stricty 
two-dimensional atomic crystals of carbon (monolayers of graphite - graphene) [21 [3] is 
promising a wealth of new phenomena and applications in technology. The observation 
of anomalous transport properties, and, most exciting, the discovery of substantial 
field effect and magnetism at room temperature allows one to envisage graphene as 
a reasonable replacement of nanotubes in electronic applications |H |5] . 

In the case of isolated graphene, the electronic states near the Fermi level can be 
described in a simple manner [6j. By symmetry, the lower and upper bands touch at the 
corners of the hexagonal Brillouin zone. In the vicinity of these points, the dispersion 
relation is isotropic and linear, and the density of states at the Fermi level is strictly 
zero. Using the tight-binding approximation for the nearest neighbour interaction in 
the honeycomb lattice, an effective long-wavelength description of electronic states in 
graphene can be written in terms of a continuum model which is based on the Dirac- 
Weyl equation for massless electrons in 2+1- dimensional space-time with the role of 
speed of light c played by Fermi velocity v ~ c/300 [3 [8]. The one-particle Hamiltonian 
operator of the model takes the form 

H = -ihv{a^d^ + a^dy), (1) 

where the 4x4 matrices can be chosen in the form p] 



a^{j = 1, 2, 3) are the Pauli matrices. The one-particle electronic wave function possesses 
4 components: 

1p = (^A+, i^A-, -ipB-f , (3) 

where subscripts A and B correspond to two sublattices of the honeycomb lattice and 
subscripts + and — correspond to two inequivalent Fermi points. The linear dispersion 
law, E = |k|, results in the density of states 

^(^) = (4) 

which is the density of states of two-dimensional ultrarelativistic Fermi gas (here S is the 
area of the graphene sample). Consequently, the thermal average charge of electronic 
excitations in graphene, 

oo 

Q(T) = -| j dET{E)te.nh(^^pEy (3 = {kBT)-\ (5) 

— oo 

is evidently zero, because Eq.(4) is even in energy [ks is the Boltzmann constant). 
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Topological defects in graphene are disclinations in the honeycomb lattice, resulting 
from the substitution of a hexagon by, say, a pentagon or a heptagon; such a disclination 
rolls up the graphitic sheet into a cone. More generally, a hexagon is substituted by a 
polygon with 6 — sides, where A^^ is an integer which is smaller than 6. Polygons 
with Nd > {Nd < 0) induce locally positive (negative) curvature, whereas the graphitic 
sheet is fiat away from the defect, as is the conical surface away from the apex. In the 
case of nanocones with A^^^ > 0, the value of A^^^ counts the number of sectors of the value 
of 7r/3 which are removed from the graphitic sheet. If A^^ < 0, then — A^^ counts the 
number of such sectors which are inserted into the graphene sheet. Certainly, polygonal 
defects with Nd > 1 and Nd < —1 are mathematical abstractions, as are cones with 
a pointlike apex. In reality, the defects are smoothed, and A^^ > counts the number 
of the pentagonal defects which are tightly clustered producing a conical shape; such 
nanocones were observed experimentally [ID]. Theory predicts also an infinite series 
of the saddle-like nanocones with — A^^ counting the number of the heptagonal defects 
which are clustered in their central regions. Saddle-like cones serve as an element which 
is necessary for joining parts of carbon nanotubes of differing radii and for creating 
schwarzite [H], a structure appearing in many forms of carbon nanofoam [12]. As 
it was shown by using molecular-dynamics simulations [13], in the case of A^^ < —4, a 
surface with a polygonal defect is more stable than a similarly shaped surface containing 
a multiple number of heptagons; a screw dislocation can be presented as the Nd ^ —oo 
limit of a 6 — A^^-gonal defect. 

In the present paper we shall consider the influence of topological defects in 
graphene on its electronic properties at finite temperatur^. The effects of the variation 
of the bond length or the mixing of vr- with cr-orbitals caused by curvature of the lattice 
surface are neglected, and our consideration, focusing on global aspects of coordination 
of carbon atoms, is based on the long-wavelength continuum model originating in the 
tight-binding approximation for the nearest neighbour interactions. We employ a general 
theory of planar relativistic fermionic systems in the singular vortex background [T9|[20]: 
its version for the case of massless fermions was elaborated in [2TI l22l [23] . 

2. Hamiltonian, resolvent and density of states on a graphitic nanocone 

The Dirac-Weyl Hamiltonian for electronic excitations on a graphitic nanocone with a 
pointlike apex takes the form 



H 



ihv [a^dr + c?r 



1 



(6) 



where 




(7) 



% The case of zero temperature was discussed earlier in works [TU [151 [13 [IZl [HI [S] • 
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is the pseudospin, r and if are the polar coordinates centred at the apex of the cone 
with metric 



2„2 



= (1 - ^) r 



and 



7] = NJQ. 



(8) 



(9) 



The electronic wave function on a graphitic nanocone obeys the Mobius-strip-type 
condition [9]: 



V'(r, v? + 27r) = - exp (-i^NdR^ ip{r, ip), 



where 



R 



a2 
-a' 



(10) 



is the operator exchanging sublattices, as well as inequivalent Fermi points, and 
commuting with Hamiltonian (6). 

By performing singular gauge transformation 



one gets the wave function obeying usual condition 

-^'{r, ip + 27t) = -ip'ir, ip), 
in the meantime, Hamiltonian (6) is transformed to 



H' = e'^He 



Q Tj -i^ _ 



ihv < a^dr + c?r 



2„-l 



(12) 
(13) 

,(14) 



where relation (9) is recalled. Thus, a topological defect in graphene is presented by a 
pseudomagnetic vortex with flux through the apex of a cone with deficit angle 

iVdVr/3. 

By performing a unitary transformation which diagonalizes i?, one gets Hamiltonian 
in the block-diagonal form 



H" 



Hi 




where 



hv < ia dr — a r 



1^-1 



-cr 
2 



(15) 



, s = ±l. (16) 



Let us consider the kernel of resolvent (Green's function) of Hamiltonian (16) 

27r 



[r, r ) awir, r 



22 V 5 



where Z is the set of integer numbers, and is a complex parameter with dimension of 
energy. Off-diagonal radial components are expressed through the diagonal ones, while 
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the latter satisfy the second-order differential equations. All radial components behave 
asymptotically at large distances as outgoing waves, i.e. as e'^^^ {2tt ^/r)~^ at r ^ cxo 
and as e*^'"'(27rv^)~^ at r' — > oo, where k = ^/oJ^ {hv)~^ and a physical sheet for the 
square root is chosen as < Arg/c < vr (ImA; > 0). In the cases of A^^^ = 3, 4, 5 all 
radial components are regular at small distances, i.e. at r — > and at r' —>■ 0. In 
all other cases there are radial components that are irregular at r = and r' = for 
some values of n. This is related to the fact that in these cases Hamiltonian (16) is not 
essentially self-adjoint and is characterized by the nonzero deficiency index. According 
to the Weyl-von Neumann theory of self-adjoint operators (see, e.g. [21]), a procedure of 
the self-adjoint extension is implemented yielding the condition for the irregular radial 
components, which depends on the set of self-adjoint extension parameters. It should 
be emphasized that only irregular radial components can produce a piece in the density 
of states, which is odd in energy. Thus, if charge (5) is nonvanishing, then it depends 
on the set of self-adjoint extension parameters. 

Let us consider the possibility of charge generation in the cases of = 
2, 1, —1, —2, —3, —6, when Hamiltonian (16) is characterized by deficiency index (1, 1), 
and there is only one self-adjoint extension parameter - 0. The irregular diagonal radial 
components are given by expressions 

a\,'"{r;r') 



X 



^ ' ^ 2(1 -r]) (sini/,^ + cosz/<^e^-^^) 
(1) 



where 



9{r — r')H_p{kr)[smh'^J_F{kr') + cos J F{kr')] + 
+e{r'-r)[smu^J^Fikr) + cosi^^JF(A;r)]i7i^|(fcr')}, (18) 

^22 [r,r ) - _ ^^.^ ^ v^e'^^) ^ 
X |6l(r - r')H^^lp{kr)[sm.iJ^Ji_F{kr') - cosu^J^i+F{kr')] + 
+e{r' - r)[smu^Ji_F{kr) - cos iy^J.i+F{kr)]H[%{kr')'j , (19) 

^ f [sgn(Ar,) - 1] N, = 2, 1, -1, -2, -3, 
' ^ -2s, AT, = -6, ^ ' 

[3-3sgn{N,) + N4{Q-Nd)-\ AT, = 2, 1, -1, -2, -3, 

1/2, N, = -6, ^''^ 

Jfj,{u) is the Bessel function of order /i, H^\u) is the first-kind Hankel function of order 
^' 

' 1, u>0 
0, M < 



parameter i^^ is related to G in the following way 

tan.. = -- — tan - + -), (22) 
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M is the scale-invariance-breaking parameter of dimension of mass. The diagonal radial 

components for n ^ n^. arc regular. 

Computing the contribution of the irregular radial components to the functional 
trace of the resolvent, we get 



[Tr{H,-uj)-\^^^ = {l-ri) J drr 



2F - 1 



all [f^f) + [r,r) 



^ 7T- (23) 

where Tr . . . = J d'^Xy^tT{x.\ . . . |x) and tr denotes the trace over pseudospinor indices 
only. Using the relation between the resolvent trace and the density of states, 

t{E) = \Tr{H -E- iO)-^ - Tr(i/ -E^ iO)"^] , (24) 

we get the expression for the piece of the density of states, which is odd in energy: 

2(2F - 1) sin(F7r) 
T^ddi-C/j = ^ X 

(a)""tan(f + f) + (ja)"^%ot(| + f) 

^ / \ 2('2F-l~) / \ 2('l— 2F1 ' V / 

{^) 'tanMf + J)-2cos(2f.)+(-ia)' 'cot^f + J) 

where the summation over s = ±1 is performed. Although quantity (25) is negligible 
as compared to the ideal gas contribution (4), lacking the factor of the sample area, it 
provides the generation of nonzero charge on some graphitic nanocones. 

3. Induced charge 

Substituting (25) into (5), we get the thermal average of charge, 

oo 

Q{T) = 2e(l - 2F)^^ / ^ tanh ( ^-uPMv'\ x 

n J u \2 J 

.^^-4an(f + f)+.^-^^cot(f + f) 
m2(2F-i) tan^ (I + I) - 2 cos(2F7r) + u^(^-^F)cot^ (f + f ) ' 
Taking into account relation (24), one can present thermal average (5) in the form 

^^^^ ^ "i/ |^Tr(i/-a;)-Hanh(i/3u;), (27) 
c 

where C is the contour consisting of two coUinear straight hues, (— oo + iO, +oo + iO) 
and (+00 — iO, — oo — iO), in the complex cj-plane. By deforming contour C to encircle 
poles of the hyperbolic tangential function on the imaginary axis, one gets 

Q(T)^~J2^(H-'''rnr\ (28) 
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where = (2m + l)7r//3. Using (23) and summing over s 



±1, we get 



X 




(2m + l)7r 



2F 



tan I - + - 



+ 



(2m + l)7r 

(3Mv^ 



2{1-F) 




Representation (26) can be regarded as the Sommerfeld- Watson transform of the infinite 
sum representation (29). 

Thus, the charge is zero for graphitic nanocones with A^^^ = 2, —2, —6, since in these 
cases one has F = 1/2, see (21). Note that one has: F = 1/5 for = 1, F = 5/7 for 
Nfi = —1, and F = 1/3 for A^^^ = —3. Therefore the charge, if any, for the one-pentagon 
{Nd = 1) defect is of opposite sign to that for the one- heptagon (A^^^ = — 1) defect and 
is of the same sign to that for the three-heptagon (A^^ = —3) defect. 

In the high-temperature hmit {(3 —>■ 0) the charge tends to zero as inverse power of 
temperature: 



:i-2 



-2+2F 



4e 

oo) = — (1 -2F) X 



)C(2 



2F) tan — + - 



X < 



(l-2--)C(2F)cot(- + - 







TT 



4^ 



TT 
2F-1 



1-2F 



TT 



0<F<- 



-<F<1 



(30) 



where ({u) is the Riemann zeta function (see, e.g. 
the charge tends to finite value [9]: 

g(0) = esgno[(l-2F) cos 9], 

where 

' 9{u)-9{-u) 



0, 



[25]). In the zero-temperature hmit 

(31) 

u = 



Thus, there is a remarkable equation relating the ground state charges which are induced 
by the one-pentagon, one-heptagon and three-heptagon defects: 

g(o)u,=i = -g(o)u,=_i = Q(o)u,=_3. (32) 

As it is shown in [9J, the induced ground state charge is accumulated around a defect. 
This suggests that the same happens for the temperature-induced charge as well, and 
the local density of states is modified in essential way in the vicinity of a defect. 

We conclude that the theory predicts the dependence of the induced charge on the 
self-adjoint extension parameter in the cases of graphitic nanocones with one pentagon, 
one heptagon and three heptagons at the apex. This is due to the fact that the density 
of electronic states attains an additional piece which is odd in energy, see (25). Actually, 
there are three possibihties: cos 9 > 0, cos 9 < 0, and cos 9 = 0. The question of which 
of the possibilities is realized has to be answered by future experimental measurements. 
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